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Peripheral Blood Monocyte Gene Expression Profile 
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Novel biomarkers of disease progression after type 1 diabetes 
onset are needed. We profiled peripheral blood (PB) monocyte 
gene expression in six healthy subjects and 16 children with type 1 
diabetes diagnosed ~3 months previously and analyzed clinical 
features from diagnosis to 1 year. Monocyte expression profiles 
clustered into two distinct subgroups, representing mild and se- 
vere deviation from healthy control subjects, along the same 
continuum. Patients with strongly divergent monocyte gene ex- 
pression had significantly higher insulin dose-adjusted HbA lc 
levels during the first year, compared with patients with mild 
deviation. The diabetes-associated expression signature identi- 
fied multiple perturbations in pathways controlling cellular me- 
tabolism and survival, including endoplasmic reticulum and 
oxidative stress (e.g., induction of HIF1A, DDIT3, DDIT4, and 
GRP78). Quantitative PCR (qPCR) of a 9-gene panel correlated 
with glycemic control in 12 additional recent-onset patients. The 
qPCR signature was also detected in PB from healthy first-degree 
relatives. A PB gene expression signature correlates with glycemic 
control in the first year after diabetes diagnosis and is present in at- 
risk subjects. These findings implicate monocyte phenotype as 
a candidate biomarker for disease progression pre- and post- 
onset and systemic stresses as contributors to innate immune 
function in type 1 diabetes. Diabetes 61:1281-1290, 2012 




The most promising stage for therapeutic inter- 
vention in type 1 diabetes is presently at diag- 
nosis, with the goal of prolonging remission and 
preserving residual (3-cell mass. However, the 
rate of type 1 diabetes progression is heterogeneous and 
unpredictable. The quality of glycemic control in early dis- 
ease has far-reaching implications for disease progression 
and the development of diabetes complications (1). Youn- 
ger age, presentation with diabetic ketoacidosis, a history of 
viral infection, or low adiponectin (a marker of insulin re- 
sistance) have been associated with a low frequency of 
remission (2-4), but no biomarkers offering insight into the 
immunological basis of disease progression have so far 
been identified. Such biomarkers would also assist trial 
design. For example, immunotherapies in recent trials 
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temporarily delayed disease progression, but treated pa- 
tients eventually deteriorated on a similar trajectory to 
control subjects (5). 

The versatile cells of the innate immune system, in- 
cluding peripheral blood (PB) monocytes and their differ- 
entiated tissue progeny, macrophages, and dendritic cells 
(DCs), are powerful sensors of stress and danger and 
mediators of tissue repair. In their role as antigen pre- 
senting cells, they orchestrate adaptive immune responses 
to both endogenous and pathogen-mediated host damage. 
Paradoxically, innate immune activation may play both 
pathological and protective roles in type 1 diabetes, de- 
pending on disease stage and context (6). At a molecular 
level, there is substantial evidence for dysregulated innate 
immunity in type 1 diabetes (7-10), both before and after 
disease onset. Activated macrophages and DCs are ob- 
served in rodent pancreatic islets before lymphocyte re- 
cruitment and are thought to be the primary source of the 
proinflammatory cytokines tumor necrosis factor (TNF), 
interleukin (IL)-6, and IL-ip that are implicated in type 1 
diabetes pathogenesis (11). Macrophages are also found in 
islets from recent-onset type 1 diabetic patients (12). In 
contrast with their proinflammatory roles, macrophages 
also have immunoregulatory and tissue repair roles, which 
may be critical in restraining (3-cell destruction. Diabetes 
was exacerbated in mice that lack colony stimulating 
factor (CSF)-l-dependent pancreatic macrophages, whereas 
macrophages promoted (3-cell survival and function in a 
CSF-l-dependent manner (13,14). 

We hypothesized that PB monocyte phenotype and 
function constitute sensitive indicators of the prediabetes 
microenvironment, making them ideal biomarkers of innate 
immune status in type 1 diabetes. Previous PB screening 
approaches identified proinflammatory signatures, although 
these studies were in mixed cell populations and were 
frequently complicated by the proinflammatory effects 
of hyperglycemia (8,9). We characterized monocyte gene 
expression in insulin-treated type 1 diabetic patients 3 
months after disease onset and analyzed clinical features 
from diagnosis to 1 year. We report a PB monocyte ex- 
pression profile, which clinically stratifies glycemic control 
within the first year after type 1 diabetes diagnosis. 

RESEARCH DESIGN AND METHODS 

Study population. Sixteen children with recently diagnosed type 1 diabetes 
(~3 months postdiagnosis; Table 1) were recruited at the Mater Children's 
Hospital during 2009 and 2010. Informed consent was obtained from all par- 
ticipants, and the studies were approved by the Mater Health Services and 
Princess Alexandra Hospital Human Research Ethics Committees. Consenting 
adult blood donors were used as healthy control subjects for the microarray 
study (33% female, mean [SD] age 32 years (8.2 years), BMI 24.5 [1.2]), whereas 
age-matched children without a family history of autoimmune disease recruited 
from operation list patients awaiting nonurgent surgery served as control subjects 
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TABLE 1 

Characteristics of group A and group B type 1 diabetic patients at diagnosis 





Group A 


Group B 


P 


n 


9 














Age at diagnosis (years ± SD) 


11.5 ± 2.9 


11.6 ± 2.8 


0.94 


olvll at diagnosis \ ± &u) 


on o -i- a i 
ZU.Z ±0.1 




U.yo 


Number of patients with >1 islet autoantibody (%) 


7 (78) 


1 f~\ C\(X\ 

( (1UUJ 


U.4o 


Mean glucose at diagnosis (mmol/L ± SD) 


24.5 ± 6.0 


26.4 ± 12.1 


0.69 


Mean %HbA lc at diagnosis (± SD) 


13.2 ± 0.9 


13.4 ± 2.3 


0.86 


Mean C-peptide at diagnosis (nmol/L ± SD) 


0.19 ± 0.1 


0.13 ± 0.03 


0.17 


Mean 25-OH-vitamin D at diagnosis (nmol/L ± SD) 


93.2 ± 37.1 


95.0 ± 25.6 


0.91 


Mean HC03" at diagnosis (mmmol/L ± SD) 


15.2 ± 8.3 


18.3 ± 6.3 


0.43 


Number of patients with celiac disease 


1 


1 


1 


Number of patients with high-risk HLADQ genotype** 


5/7* 


3/5* 


1 



Biometrics and clinical parameters recorded at diagnosis of type 1 diabetes. P values calculated using two-tailed unpaired t test (discrete 
variables) or Fisher exact probability test (categorical variables). HC03~, bicarbonate; DKA, diabetic ketoacidosis. *No genotyping data for 1 
of 9 group A and 2 of 7 group B patients. **DQB1*0201, DQBP0302 and/or DQAP0501, DQAP301. 



for flow cytometry and peripheral blood mononuclear cell (PBMC) quantitative 
PCR (qPCR) (mean age [SD] 10 years [4.2 years]). Blood samples were obtained 
from first-degree relatives (FDR) of type 1 diabetic patients with informed 
consent (mean age [SD] 9 years (3.9 years). Insulin dose-adjusted glycosylated 
hemoglobin (HbAi c ) (IDAAlc) was calculated as %HbA lc + 4 X (insulin dose 
[units/kg/day]) (15). Diabetic ketoacidosis (DKA) at presentation was defined 
as <22 mmol/L plasma bicarbonate. 

Monocyte purification, expression profiling, and qPCR. CD14 + mono- 
cytes were purified from Ficoll-purified PBMC by positive selection (MACS, 
Miltenyi). Total RNA was purified from CD14 + monocytes using RNeasy Mini 
columns, with on-column DNase digestion (Qiagen). RNA from the donors 
selected for microarray (500 ng; Agilent Bioanalyzer RNA Integrity Score >8) 
was amplified using the Totalprep RNA amplification kit (Ambion). Amplified 
RNA was hybridized to Illumina HT-12 v3.0 whole genome beadarrays, 
according to the manufacturer's instructions. Beadarrays were scanned on an 
Illumina iScan, and raw data were exported from BeadStudio (Illumina) for 
further analysis. Gene expression data have been deposited in the gene ex- 
pression omnibus (GEO) database (GSE33440, Reviewer Access: http://www. 
ncbi.nlm.nih.gov/geo/query/acc.cgi?token=dvezxukmesmsohe&acc=GSE33440). 
Data were variance stabilized and quantile normalized using the Lumi package 
in R (Bioconductor, http://www.bioconductor.Org/packages/2.0/bioc/html/lumi. 
html), followed by clustering and differential expression analysis in GenespringGX 
(Agilent). Probes that did not achieve an Illumina Detection P value = 0 in >1 
sample were omitted from the analysis. Genes whose expression was signifi- 
cantly correlated with plasma glucose or %HbA lc were identified by quantitative 
trait analysis in BRB-ArrayTools, developed by Dr. Richard Simon and the 
BRB-ArrayTools Development Team (Pearson correlation, P < 0.001). Probes 
identified as significantly up- or downregulated in T1D-B monocytes were 
analyzed for ontology enrichment using Database for Annotation, Visualization 
and Integrated Discovery (DAVID; http://david.abcc.ncifcrf.gov/). Functional 
Annotation Charts generated using Level 4 Bioprocess Gene Ontology terms, 
with the Illumina HT-12 genome as background, were manually curated to remove 
redundant annotations. OligodT-primed cDNA was synthesized from 1 (xg RNA 
(Quantace). qPCR was performed using Sensimix SYBR mastermix (Quantace) 
on an Agilent HT-7900. Gene expression relative to hypoxanthine-guanine 
phosphoribosyltransferase was calculated using the ACt method. Median- 
normalized qPCR data were clustered in GenespringGX. qPCR primer sequences 
used in this study appear in Supplementary Table 1. 

Flow cytometry. Whole blood (50 |xL) was stained with anti-CD 14-FTTC, 
anti-CD 16-PacificBlue, anti-CX3CRl-Alexa647, or corresponding isotype 
controls (Biolegend or Becton Dickinson) and Aqua live/dead discriminator 
(Invitrogen) for 20 min. Erythrocytes were lysed with the addition of 450 |xL 
FACS Lysing Solution (Becton Dickinson) for 5 min, followed by sample fix- 
ation with 100 fxL of 1% formalin immediately before analysis on a Gallios flow 
cytometer (Beckman Coulter). Live, CD14 + cells were analyzed for CD 14, 
CD 16, and CX3CR1 expression, with gates set according to negative controls, 
which contained the isotype control for each antibody, as well as the rest of 
the antibodies in the staining panel. Kaluza (Beckman Coulter) was used for 
flow cytometry compensation and analysis. 

Statistical analysis. Except for microarray analysis, all statistical tests were 
calculated in GraphPad Prism. Comparisons between groups were made with 
unpaired t test or ANOVA; F test was used to check whether variances were 



equal, and appropriate tests were used: unpaired t test if variances were equal, 
with Welch's correction if variances were unequal. Categorical variables were 
compared using the Fisher exact probability test. 



RESULTS 

Stratification of type 1 diabetic patients by monocyte 
expression profile at diagnosis. Purified PB CD14 + 
monocytes from a total of 16 recent-onset type 1 diabetic 
patients and six healthy control subjects were expression 
profiled using whole genome microarrays (Illumina 
HT-12), in two separate, independent experiments. Because 
of the volume of blood required to study purified monocytes, 
the healthy control subjects comprised young adults. Blood 
samples for microarray analysis were taken from insulin- 
treated patients ~3 months after type 1 diabetes diagnosis, 
to avoid the acute metabolic disturbance associated with 
disease onset. A strikingly similar pattern was observed in 
both experiments, whereby a subset of monocyte expres- 
sion profiles clustered apart from the remaining patients 
and healthy control subjects (group B, Fig. 1A and E). Be- 
cause of this stark segregation, group B profiles were ana- 
lyzed separately to identify differentially expressed genes. 
In group B monocytes compared with healthy control 
subjects in both microarrays (Fig. 1C and Supplementary 
Table 2), 1,107 probes (1,015 genes) were at least 1.5-fold 
differentially expressed (Benjamini Hochberg-corrected 
P value <0.05). Between the entire diabetes cohort and 
healthy control subjects in both microarrays (Fig. ID), 320 
probes were differentially expressed. Only 38 probes were 
reproducibly differentially expressed between group A 
patients and healthy control subjects. Hierarchical clustering 
of group B-regulated probes revealed that the two putative 
type 1 diabetes subgroups represent a spectrum of devia- 
tion from healthy control subjects, along the same contin- 
uum (Fig. 2A and E). 

Extreme divergence from healthy monocyte gene 
expression correlates with early type 1 diabetes 
progression. At diagnosis, group A and group B patients 
were clinically indistinguishable in terms of HLA genotype, 
family history, clinical presentation at diagnosis (including 
symptom duration, ketoacidosis, blood glucose, 25-hydroxy- 
vitamin D, diabetes-associated autoantibodies, or the pres- 
ence of celiac or thyroid disease; Table 1 and data not 
shown). Group B patients had higher levels of HbA lc 3 and 6 
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FIG. 1. Unsupervised clustering of monocyte gene expression profiles 
identifies a subgroup of recent-onset type 1 diabetic patients. Purified 
PB monocytes from patients (D) or healthy control subjects (HC) (H) 
were expression-profiled by Illumina whole genome microarray, in two 
independent experiments (A, 1?). For each microarray, all genes detected 
in ^1 sample were clustered using the Euclidean Distance correlation. 
Monocyte expression profiles from all patients, or the subset of patients 
designated "B", were used to identify genes that were significantly 
differentially expressed (^1.5 -fold differential expression, Benjamini 
Hochberg corrected P < 0.05) between patients and control subjects 
(C,D). 



months after diagnosis (Fig. 3A) and required significantly 
higher insulin doses to maintain glycemic control in the first 
12 months after diagnosis (Fig. 3B, P = 0.0004, 2-way ANOVA). 
Approximately 30% of children diagnosed with type 1 diabetes 
experience a partial remission phase (honeymoon), during 
which their HbA lc levels and insulin requirements (used as 
a proxy for endogenous insulin production from residual 
P-cells) (15) are relatively low. An insulin dose-adjusted 
HbA lc metric (IDAAlc) <9 correlated with the definition of 
C-peptide responders in the Diabetes Control and Compli- 
cations Trial (16) and was proposed as a definition for partial 
remission. Group A and group B IDAAlc diverged around 
9 at 3-6 months postdiagnosis, and group B patients 
exhibited significantly higher IDAAlc during the first year 
of disease (Fig. 3C, P < 0.0001, 2-way ANOVA). 

To examine whether the monocyte gene expression 
profiles observed were directly related to hyperglycemia, 
we identified probes whose expression was significantly 
correlated (Pearson P < 0.001) with plasma glucose con- 
centration (n = 25), or %HbA ic (n = 70), at the time blood 
was sampled for profiling. No glucose-correlated probes 
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FIG. 2. Extremely divergent monocyte expression profile in a snbgronp of 
recent-onset type 1 diabetic patients. Genes that were differentially 
expressed between group B monocytes and healthy control subjects (HC) 
(ji = 1,1^7) in two independent microarrays were clustered using the Eu- 
clidean Distance correlation, as were the experimental samples (A, 1?). 



and 9 of 70 HbAi c -correlated probes were differentially 
expressed between group B and healthy control subject 
monocytes (Supplementary Fig. 1). These analyses suggest 
the innate immune status of group B monocytes is not solely 
driven by hyperglycemia. 

The group B monocyte profile is a signature of 
disturbed metabolic homeostasis. The correlation be- 
tween monocyte gene expression close to type 1 diabetes 
onset and early clinical trajectory may imply a role for 
monocyte function in disease progression. Alternatively, 
or additionally, because PB monocytes are professional 
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FIG. 3. Diabetes progression in recent-onset type 1 diabetic patient 
subgroups. Plasma levels of %HbA lc (A) and total daily insulin dose 
(units/kg) (1?) were recorded at 3 monthly intervals following diabetes 
diagnosis, from which IDAAlc (C) was calculated. Mean %HbA lc , daily 
insulin dose, and IDAAlc (± SD) are shown for patients subgrouped 
according to monocyte expression profile (group A, #; group B, ■). 
Group sizes: HbA lc data: n = 8, 7, 7, and 7 for group A and n = 7, 7, 6, and 
7 for group B at 3, 6, 9, and 12 months. TDD data: n = 8, 7, 6, and 4 for 
group A and n = 6, 6, 5, and 4 for group B at 3, 6, 9, and 12 months. 
Missing data arose when patients missed a clinic visit or commenced 
a basal/bolus or insulin pump regimen. 



stress-sensing cells that circulate throughout the body, 
they may provide a window on systemic stresses, i.e., 
reflecting islet damage. We analyzed the set of 1,107 
probes that were >1.5X up- or downregulated in group B 
patients compared with healthy control subjects in both 
arrays for significantly enriched gene ontology terms. 
Terms were manually curated to remove redundant cat- 
egories comprising overlapping gene sets. Overexpressed 
genes were enriched for cellular stress and activation 
pathway genes, especially oxidative and endoplasmic re- 
ticulum (ER) stress (Fig. 4). We did not observe upregula- 
tion of classical inflammatory mediators, with the exception 
of tumor necrosis factor (TNF), interleukin (IL)-10, and 
PTGES2 (which encodes the enzyme catalyzing prosta- 
glandin E2 production), and the prostaglandin E2 receptor, 



PTGER (Supplementary Table 2 and Fig. 5), whereas mul- 
tiple genes involved in regulation of immune responses 
were downregulated in group B monocytes (Fig. 4). These 
included sensors such as TLR4, TLR5, IFI16, NLRX1, and 
NLRC4/IPAF and effectors such as CASP1 (Supplementary 
Table 2). Ten genes in key pathways differentially ex- 
pressed by group B monocytes were selected for qPCR 
validation. Their expression was quantified in eight patients 
profiled previously by microarray and 12 additional recent- 
onset type 1 diabetic patients. Hierarchical clustering of 9 
robustly detected genes segregated group B patients from 
group A and healthy control subjects, whereas unclassified 
patients' monocytes (D) revealed a spectrum of gene ex- 
pression, similar to the microarray analysis (Fig. 5A). IDAAlc 
tended to increase with increasing distance from the healthy 
control subjects expression profile (Fig. 5B; 3 or 6 month 
linear regression: slope = -0.225/-0.176, P = 0.017/0.091). 
Two genes selected for validation, DDIT4 and adrenome- 
dullin (ADM), were not incorporated in the cluster analysis 
since they were only detected in monocytes from a subset of 
type 1 diabetic patients, including group B (Fig. 5C and D). 

The group B monocyte expression profile strongly sug- 
gested stress-induced metabolic disturbance, with upre- 
gulation of ER and oxidative stress responsive genes, and 
downregulation of genes involved in cellular energy and 
metabolite production, especially mitochondrial oxidative 
phosphorylation (e.g., PDHB, MDH1, IDH1, SDHC, ACLY) 
(Fig. 5 and Supplementary Table 1). Mitochondrion was 
the most significantly enriched cellular component ontol- 
ogy among group B downregulated genes (P = 2.6e-05). It 
is well established that activation of myeloid and other 
immune cells is accompanied by a switch in glucose me- 
tabolism and cellular energy production, from oxidative 
phosphorylation to anaerobic fermentation through glycol- 
ysis (17). PFKFB3, which encodes the rate-limiting enzyme 
in the glycolytic pathway, was upregulated in group B 
monocytes (Fig. 5), but glycolytic pathway genes were not 
generally induced. PDHB and PFKFB3 expression was 
strongly inversely correlated in the microarray (r = —0.86, 
P < 0.0001), which may imply a reduction in oxidative 
phosphorylation and increased glycolysis (17). In addition 
to perturbations in core metabolic pathways, there was 
also evidence of cellular activation through stress, in- 
cluding the unfolded protein response that is a conse- 
quence of ER stress (e.g., IRE1, GRP78, DDIT3, XBP1) 
(Fig. 5 and Supplementary Table 2). HIF1A, a key mediator 
of oxidative stress, was upregulated in group B monocytes, 
as were multiple HIF1A target genes, including DDIT4, 
PFKFB3, and ADM (17,18) (Fig. 5 and Supplementary 
Table 2). Thus, the striking group B monocyte expression 
profile is a mixed signature consistent with the existence of 
multiple pathological processes and/or adaptive responses 
to stress and inflammation. 

Monocyte subset balance is altered in PB of type 1 
diabetic patients and a subset of FDR. The mRNA 
encoding CX3CR1, a marker of CD16 + monocytes (19), 
was underexpressed in type 1 diabetic patients compared 
with healthy control subjects (Fig. 5), suggesting this dys- 
regulation could derive from a perturbation in PB monocyte 
subset balance. PB monocytes from children with recent- 
onset or established type 1 diabetic and healthy children 
were phenotyped in whole blood according to surface ex- 
pression of CD 14 and CD 16 (see Supplementary Fig. 2 for 
representative flow cytometric profiles). We observed a re- 
duction in the total proportion of CD16 + monocytes in PB 
of type 1 diabetic patients (not shown), which was as a 
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FIG. 4. Biological processes associated with group B differentially expressed genes. Gray bars indicate ontology enrichment within upregulated 
genes; white bars indicate ontology enrichment within downregulated genes. Associated P values are indicated by the line. NF-kB, nuclear 
factor- kB. 



result of a selective reduction in the proportion of CD14Hi/ 
CD16 + monocytes (Fig. 6A). The proportions of CD14Low/ 
CD16 + and CD14Hi/CD16~ populations were unaltered 
(Fig. 6B and C), whereas the proportion of CD14Low/ 
CD16~ cells was significantly increased in recent-onset 
type 1 diabetic patients (Fig. 6D). We observed a signifi- 
cant reduction of CD14Hi/CD16 + monocytes and an in- 
crease in CD14Low/CD16~ monocytes in FDR of children 
with type 1 diabetes, relative to healthy children (Fig. 
6E), suggesting that this finding is not associated with 
hyperglycemia. This phenotype may be monocyte-intrinsic, 
or may arise because of external sources of stress, which 
could affect p-cell reserve and PB monocyte subsets in both 
recent-onset patients and in FDR. The proportion of CD16 + 
monocytes in PB positively correlated with monocyte sur- 
face CX3CR1 expression by all CD14 + monocytes (Fig. 6F, 
Pearson's r = 0.4648, P < 0.0001), whereas CX3CR1 surface 
expression by individual monocyte populations was not al- 
tered (data not shown). Thus, the distinctive monocyte 
signature in recent-onset diabetes, which includes CX3CR1 
repression, may partially derive from a relative deficiency in 
CD14Hi/CD16 + monocytes and enrichment in CD14Low/ 
CD16~ monocytes. 



Subtle changes in cellular composition could underpin 
differential gene expression in a mixed cell population. We 
thus investigated whether group B-regulated genes typi- 
cally exhibited monocyte subset-specific expression pat- 
terns. In a meta-analysis of all 16 monocyte microarrays, 
98 genes were significantly correlated with CX3CR1 ex- 
pression (Pearson, P < 0.001). Of the 1,109 genes that 
were significantly differentially expressed in group B mono- 
cytes (>1.5X, P < 0.05), 64 were also correlated with 
CX3CR1 expression (Supplementary Fig. 3). Given the 
strong correlation between total monocyte CX3CR1 surface 
expression and %CD16 + monocytes in PB, these represent 
candidate subset-specific genes. We also analyzed two 
public microarray datasets that compared CD14Hi/CD16~, 
CD14Hi/CD16 + , and CD14Low/CD16 + monocyte gene 
expression (20,21). There was negligible overlap between 
the group B monocyte signature and the genes that were 
differentially expressed between monocyte subsets (data 
not shown). 

ER and oxidative stress markers are upregulated in 
FDR PBMC. To investigate whether elements of the group 
B monocyte signature are evident in children at risk of 
developing type 1 diabetes, we assayed gene expression in 
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FIG. 5. Quantitative PCR (qPCR) validation of monocyte gene expression profiles in recent-onset type 1 diabetes. The expression of 9 group B- 
regulated genes was quantified by real time qPCR in CD14 + monocytes, relative to hypoxanthine-guanine phosphoribosyltransferase expression. 
Relative expression values for each gene were normalized to the gene's median expression across all samples. Genes and samples were clustered 
by Spearman correlation (A). Group A and group B patients were classified previously by microarray. D denotes recent-onset patients who were 
followed from diagnosis with monocyte sampling at 3 months, but not profiled by microarray. IDAAlc corresponding to each patient profile in A 
recorded 3 (circles) and 6 (squares) months postdiagnosis (1?) is shown. Missing data arose when patients missed a clinic visit or commenced 
a basal/bolus or insulin pump regimen. Relative expression of DDIT4 (C) and ADM (Z>) in the patient cohort depicted in A and B is shown. 
Transcript detection rates appear in parentheses. HC, healthy control subjects. 



PBMC from FDR, including 14 of 43 with at least one islet 
autoantibody, compared with healthy children. Purified 
monocytes were not studied since blood sample volumes 
were limiting. Broadly expressed genes such as CX3CR1 
and HIF1A were not differentially expressed between FDR 
and healthy control subjects (data not shown). By con- 
trast, increased expression of the inducible ER stress 
(DDIT3, GRP78) and oxidative stress genes (DDIT4, 
ADM), and the proinflammatory cytokine TNF, was evi- 
dent in a subset of FDR (Fig. 1A-E). Expression levels of 
several of these genes were highly correlated to each other 
(P < 0.0001, Fig. 7F-H and data not shown), but their 
expression did not differ between autoantibody positive 
and negative FDR (not shown). Stress marker expression 
did not correlate with monocyte subset distribution (Fig. 
71 and J and data not shown), confirming our conclusion 
from the meta-analysis that the stress signature is co- 
incident with, but distinct from, altered monocyte subset 
distribution in type 1 diabetic patients and FDR. 

DISCUSSION 

There is considerable interest in developing therapies to 
protect and regenerate residual (3-cell function, to arrest 
or slow the course of disease, and to mitigate the de- 
velopment of type 1 diabetes complications (5). This 
should be achievable, since 80% juvenile diabetic patients 
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had significant residual (3-cell function during the first year 
after diagnosis (22). However, disappointing results in 
several recent Phase III immunotherapy trials in recent- 
onset type 1 diabetes have invited speculation that factors, 
such as (3-cell stress, contribute to continuing (3-cell de- 
struction after disease onset and may require treatments 
other than immune intervention (23). 

We report here that PB of children with type 1 diabetes 
is characterized by a relative deficiency of CD14Hi/CD16 + 
monocytes, an enrichment in CD14Low/CD16~ mono- 
cytes, and a striking, coordinated monocyte gene expres- 
sion reflecting altered metabolism associated with stress. 
Recent-onset patients could be reproducibly stratified on 
the basis of their monocyte gene expression profile. Severe 
deviation from a healthy monocyte expression profile, 
which we observed in —50% of patients, independently 
identified a clinical subset with reduced glycemic con- 
trol, as evidenced by higher %HbA ic and higher insulin 
requirement. We hypothesize that the two ends of the 
spectrum at type 1 diabetes onset represent different dis- 
ease trajectories, whose differences in glycemic control 
may affect future development of diabetes complications 
(1). Identification of individuals on a severe trajectory may 
provide opportunities for better clinical management or 
targeted therapies for these patients. The reason for the 
early clinical divergence in these patients is not clear, but 
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FIG. 6. Peripheral blood monocyte subset distribution in type 1 diabetic patients and their at-risk relatives. The percentage of monocytes in each of 
the 4 monocyte subsets was determined by flow cytometry analysis of surface CD14 and CD16 expression in whole PB. A-D: PB monocyte subset 
distribution was determined for healthy control subjects (HC) (circles), long-standing type 1 diabetic patients (D, squares, >12 months after 
diagnosis), recent-onset type 1 diabetic patients within 3 months of diagnosis (D[Dx], triangles). PB monocyte subset distribution in healthy 
control subjects (filled shapes) and FDR (unfilled shapes) (i?) is shown. Correlation between %CD14 + /CD16 + monocytes and monocyte CX3CR1 
surface expression (F) is shown. ***P < 0.0001; **P < 0.01; *P < 0.05. 



could be explained by two testable hypotheses: 7) patho- 
logical processes are intrinsically more aggressive, or 
regulation less effective, in group B patients than group A 
patients or 2) group A patients possess the greater (3-cell 
reserve at diagnosis, and insulin treatment enables at least 
transient restoration of function in these cells. 

Type 1 diabetes progression is most commonly assessed 
using metabolic parameters, including %HbA lc , insulin re- 
quirement, and C-peptide levels. C-peptide is the most 
common surrogate marker for declining (3-cell function, 
and thus rate of disease progression. Higher %HbA lc and 
insulin requirements were associated with rapid disease 
progression in several prospective studies, along with risk 
factors such as younger age, presentation with ketoacidosis, 
and GAD autoantibody status (4,24,25). Percent HbA ic 
levels also predict the development of microvascular com- 
plications of type 1 diabetes (26). Type 1 diabetes pro- 
gression and remission have rarely been studied from an 
immunological perspective, although increased numbers 
of regulatory T cells were associated with remission in 
one small cohort study (27). It remains to be determined 
whether monocyte phenotype represents an independent 
biomarker for disease progression or whether monocytes 
are simply stress-prone cells that reflect pancreatic dys- 
function. Notably, the mRNA encoding ADM, a vasoactive 
hormone involved in protection against oxidative stress 
associated with diabetes complications and dysglycemia 
(25,28), was upregulated in group B patients and FDR. 



Reactive oxygen species (ROS), the causative agents of 
oxidative stress, are produced in response to environ- 
mental stimuli and as by-products of cellular metabolism. 
The distinct oxidative and ER stress signature we ob- 
served in group B patients suggests these patients' mono- 
cytes are intrinsically vulnerable to stress, or that they 
exist in a stressful environment. Stress pathway activation 
has implications for monocyte differentiation and function. 
In addition to the elevated expression of oxidative stress 
responsive genes, such as HIF1A and HIF1A targets, mul- 
tiple genes involved in cellular antioxidant pathways were 
downregulated in group B monocytes (e.g., CAT, G6PD, 
OXR1, PRDX1, and PRDX3), suggesting protective mech- 
anisms may be impaired. Oxidative stress is strongly im- 
plicated in type 1 diabetes progression, particularly in the 
development of complications, in both humans and animal 
models (29,30). Oxidative stress is directly associated with 
(3-cell toxicity and decreased insulin secretion, since 
P-cells have intrinsically low antioxidant defenses (29), but 
systemic oxidative stress may also play a role in promoting 
hyperglycemia and diabetes progression. Monocytes and 
neutrophils from type 1 diabetic subjects generated more 
ROS (31), and macrophage-derived ROS were found to be 
critical for the induction of autoreactive T cells and di- 
abetes development in NOD mice (32). 

Oxidative and ER stress are closely linked. Correct ER- 
dependent protein folding is essential for the localization 
and function of secreted and membrane-bound proteins. 
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FIG. 7. Elements of the type 1 diabetes group B expression signature detected in healthy FDR PBMC. Expression of DDIT3 (A), GRP78 CB), DDIT4 
(C), ADM (2>), and TNF (2?) in PBMCs from type 1 diabetic FDR or healthy children (HC) was quantitated by PCR. Correlation between DDIT3 
and GRP78 expression levels (F) is shown. Correlation between DDIT4 and ADM expression levels (G) is shown. Correlation between TNF and 
ADM expression levels (7T) is shown. Correlation between DDITS expression and %CD14Hi/CD16 + (J) and %CD14Low/CD16" (J) monocytes in PB 
is shown. P values for group comparisons (A-D) were generated by Student t test (with Welch's correction where variances were significantly 
different, as determined by F test). Pearson's correlation coefficient (r) and associated P values are shown in E-G. 
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The UPR is a transcriptional program activated in re- 
sponse to a variety of stressors that disturb ER homeo- 
stasis, including ROS, infection, lipotoxicity, and impaired 
calcium homeostasis (33). Although first characterized as 
a physiological response regulating cell survival under 
stressful conditions, emerging data demonstrate extensive, 
bidirectional cross-talk between the UPR, inflammatory, 
and metabolic pathways (33). For example, ER stress re- 
duced proinflammatory responses in monocytes from type 
2 diabetic patients (34) and nonimmune cells (35). By 
contrast, ER stress has been shown to synergise with in- 
flammatory stimuli in macrophages, specifically enhancing 
IL-23, IL-8, and IL-6 production, but not affecting IL-12 or 
IL-1(3 (36,37). The duration and context of ER stress may 
thus influence the balance of monocyte inflammatory 
mediators. Notably, TNF was upregulated in group B 
monocytes, whereas other inflammatory mediators and 
signaling molecules were downregulated, including several 
components of the IL-ip processing pathway that have 
previously been positively associated with type 1 and type 
2 diabetes (38). 

Persistent ER stress leads to apoptosis and inflam- 
mation, in the event the UPR fails to restore ER homeo- 
stasis. The induction of DDIT3, a key UPR proapoptotic 
effector, in group B monocytes suggests extreme or pro- 
longed ER stress. Because mitochondria are necessary for 
stress adaptation in the early phase of the UPR, mito- 
chondrial defects should render cells more vulnerable to 
ER stress-induced death (39). The striking downregulation 
of mitochondrial component and pathway genes in group 
B monocytes may affect susceptibility to stress, immune 
function (40,41), and cellular energy metabolism (42). 
Genes regulating apoptosis were also enriched in group B 
monocytes, which may be evidence of a stress-induced 
challenge to survival in these cells. The detection of stress 
elements of the group B monocyte signature in unfractio- 
nated PBMC from FDR confirms these gene expression 
changes are not simply a result of hyperglycemia and raises 
the possibility that oxidative and ER stress also contribute 
to disease development. Future studies will determine 
whether the PBMC stress gene expression signature derives 
solely from monocytes or is common to multiple cell types 
and will clarify the relationship between monocyte pertur- 
bation, diabetes risk status, and disease development. 

Monocytes and their differentiated progeny, macro- 
phages and DCs, play both proinflammatory and anti- 
inflammatory, tissue protective roles. Type 1 diabetic 
patients and their FDR exhibited a reduced proportion 
of CD14Hi/CD16 + intermediate monocytes in PB and an 
apparently concomitant increase in the CD14Low/CD16~ 
monocyte subset. Intriguingly, perturbations in PB mono- 
cyte subsets did not correlate with the induction of oxi- 
dative and ER stress genes in FDR PBMC, suggesting the 
two phenotypes are independent. Expansion of the CD16 + 
population has been observed in many inflammatory dis- 
eases, from infection to autoimmunity (19). A decrease in 
CD14Hi/CD16 + monocytes was associated with cardio- 
vascular events in patients with kidney disease (43,44) and 
with asthma severity (45). By contrast, increased CD14Hi/ 
CD16 + monocyte counts correlated with good clinical 
outcomes after stroke, potentially because of their roles in 
tissue repair and angiogenesis (46). CD16 + monocytes 
have also been implicated in promoting immune tolerance 
(47,48). These two key monocyte functions may be 
mechanistically linked in type 1 diabetes, since defective 
clearance of dying cells is also hypothesized to promote 



inflammation and autoimmunity (49). Wound healing is 
defective in humans and mice with diabetes, due in part to 
an inability to efficiently clear apoptotic cells (50). Thus, 
a relative deficiency in CD16 + monocytes may compromise 
tissue repair and/or autoimmune regulation at sites of in- 
flammation, including the pancreas, which would likely 
promote disease progression and complications. 

In conclusion, we demonstrate that monocyte gene ex- 
pression distinguishes type 1 diabetic patients at disease 
onset. This stratification correlates with early disease tra- 
jectory, which we hypothesize has implications for the rate 
of disease progression and the development of complica- 
tions. The coordinated expression signature in rapidly 
progressing type 1 diabetes suggests an underlying high- 
level susceptibility or exposure to physiological stress. 
Further studies will clarify the impact of the group B 
monocytes phenotype on antigen-presenting cell function, 
and thus T cell responses in autoimmune diabetes. 
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